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Abstract 

The intermittent dynamics of the turbulent GOY shcU-modcl is char- 
acterised by a single type of burst-like structure, which moves through the 
shells like a front. This temporal structure is described by the dynam- 
ics of the instantaneous configuration of the shell-amplitudes revealing a 
approximative chaotic attractor of the dynamics. 



1 Introduction 

One of the main goals in current turbulence research is to understand the ef- 
fect of intermittency in turbulence. It has long been known that intermittency 
produces corrections to the classical Kolmogorov —5/3 scaling law and to other 
moments of the energy spectrum in the inertial range ^. Still very little is 
known about the intense intermittent structures found in turbulent flows 

Over the last decade turbulent shell-models have been studied intensively 
because of their simplicity and excellent agreement of their statistics in com- 
parison with that of experimentally measured turbulence |Q, ^, ^ |7|, |8|, |9| . In a 
way these models reproduce the statistics far better than expected from their 
simplicity, so the general idea has been to reveal the nature of the dynamics of 
these models and then afterwards relate this experiences to the full problem of 
turbulence. 

While the main approach to shell models has been statistical, much can 
be learned from the study of the temporal structures naturally arising in the 
model. Because of the intense dynamics during these temporal events they are 
called bursts. Only recently have these temporal structures been thoroughly 
studied together with the nature of their creation |ll|. Among the large 
numbers and types of different shell-models the present work is based on the 
successful GOY shell- model [|, |] . 

The paper is organised as follows: The first part gives a detailed descrip- 
tion of the bursts of the standard GOY-model. The second part shows that an 
approximative chaotic attractor of the model exists expressed by the collective 
dynamics of the neighbouring shells. 
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2 The GOY-model 

All shell-models simulate the flow of energy through wavenumber space in fully 
developed turbulence. The models consist of a system of coupled ordinary 
differential equations where the energy is injected into the low wave-numbers 
by a constant forcing term. The energy then cascades up to the high wave- 
numbers by means of a coupling term where it is dissipated away by a viscosity 
term. 

2.1 Construction 

In the GOY model wave-number space is divided into N separate shells with 
characteristic wavenumbers kn = koX"" (A = 2) where n = 1, . . . , N and ko is a 
constant determining the smallest wavenumber in the model. Each shell is as- 
signed a complex amplitude which can be imagined as the velocity difference 
on a scale in = ^/kn- By assuming conservation of phase space, energy and 
helicity and interactions among the nearest and next nearest neighbour shells, 
one can arrive at the following set of evolution equations |^ |9| 

+ I'k^ Un = ikn{Un+lUn+2 ~ '^''^n-l'^n+l ^ ^n-2^n-l) + f^n,4: (1) 

with boundary conditions m_i = uq = un+i = un+2 = 0, and constant forcing 
/ on the fourth shell. 

The set (|^) of N coupled ordinary differential equations is numerically in- 
tegrated by standard techniques [0]. In the simulations, we use the following 
standard values: 5 = 1/2, N = 19, u = 10-^ k^ = 2"^, / = (1 + i) * 0.005 as 
found in earlier work [^, ^, 0, |9| . 

2.2 Conservation laws, fixed-points and invariance 

The strength of the shell-models relates to their simplicity in construction, 
which can be justified by imposing the same conservation-laws and invariants 
as the Naiver-Stokes Eq. i.e. the basic equations governing fluid dynamics. As 
for other shell-models the GOY-model exhibits these conservation laws in the 
absence of forcing and viscosity (/ = = 0) reducing the right side of Eq.|^ to 
the coupling-term. 

The conservation of phase-space is enforced as the coupling term does not 
contain Un- The remaining conserved quantities are quadratic i.e. they can be 
written in the form: = X^^nl^nP- Using the relation \{u\) = UnUn and 
inserting Qq, in the model, the coupling term becomes three terms of three 
successive amplitudes multiplied together. Comparing these three terms gives 
the following relation on a: 1 — (52" — (1 — (5)2^° = with two solutions: a = 
and a = — \n\{5 — 1). The first (a = 0) corresponds to the conservation of 
energy while the other solution corresponds to helicity conservation in the case 
of three dimensional turbulence {5 < \) and to enstrophy conservation in the 
case of two dimensional turbulence {5 > 1) 0, ^ 

The fact that the coupling term multiplied by u„ gives three terms all similar 
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within pre-factors and a displacement in indices makes the dynamics of the 
model invariant to the following changes in the complex phase: 



Vj'j-i > 6 Ufi^ 

Un+i — > e*^~"ti„+i, where n modulus 3 = 1 (2) 

Un+2 e*"'^n„+2, 

and a and /3 are free parameters. This invariance affects not only the phases 
but also the dynamics of the model since every third shell tend to follow the 
same behaviour. 

Thinking of the GOY-model as a dynamical system a basic thing to study is 
the fixed-points of the model: iin = 0, n = 1 . . . N . Requiring again an inviscid 
and unforced model (/ = = 0) gives two non-trivial scaling fixed-points: 
Un = g{n) with = | and 2 = ^(1 — \n.\{5 — 1)) where g{n) is an arbitrary 



function of period three in n coming from the invariance of the model. The 
first fixed-point: Un ~ kn^^^ g{n) corresponds to the Kolmogorov — | scaling- 
law and will be called the Kolmogorov fixed-point while the other solution result 



in an alternation of the amplitudes 11]. In spite the simplicity of the these 



fixed-points they play a crucial role in the later analysis of the model. 



3 Dynamics of the model 

At large time-scales the dynamics of the model may seem stochastic but as 
the time-span refines distinct spikes emerge and in the end the dynamics is 
noiseless and fully resolved even during the most dramatic changes. To observe 
the general behaviour we monitor |Mn(i)| as a function of time and because of 
large variations in magnitude it is shown in a semi-logarithmic plot in Fig.|l[ 
The higher shells has the lowest absolute value and fastest variations while the 
lower shells vary over a longer time-span. 
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Figure 1: A typical evolution of the norm of the amplitudes |mi| . . . |ui8| 

Two main features strike from Fig|l|: All the higher shells evolve in a syn- 
chronised manner and the evolution follows a general pattern of strong bursts 
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interchanged by oscillatory relaxation. Bursts occur randomly in time with 
great variations in their strength. 

3.1 Organisation of shell-dynamics 

The synchronisation of the different shells is a result of the coupling between the 
shells making the model self-organise into the types of behaviour seen in Fig.|l[ 
The organisation in the model is showed most clearly by the local two-point 
correlation function measuring how the dynamics of a given shell is correlated 
to its neighbourhood of both shells and in time. It is defined using the following 
two shortcuts 

= Unoit), Ua = C/„o+An(i + At) 

by: 

r(At, An) = C{Uo, Ua) = - 

and where the averages are taken over time. 

The information gained from the |r(At, An)| is divided into two parts: First 
only the norm of the complex amplitudes is correlated replacing Uq and C/a by 
their norms. This is shown in the left part of Fig.Q as a contour-plot with dark 
as the strongest normalised correlation. Second the full complex amplitudes 
are correlated and showed in the same manner in the right part of Fig.^. Both 
correlations have no = 15 and averaged over 40.000 n.u. which correspond 
roughly to a time-span of approximately 4000 successive bursts. 




Figure 2: The two-point correlation in shells and time first for the norm of the 
amplitudes and second for the pure amplitudes, based at the 15'th shell. 

The left plot show that all the amplitudes in the model is strongly correlated 
from the forcing at the 4 shell up to the highest shells. This strong correlation 
is due to the organisation of the amplitude dynamics during both bursts and 
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the succeeding strong oscillations. The same plot also shows the motion of the 
burst through the shells by the time-shift of the correlation-peeks for increasing 
An. When taking the amplitude phases into account the correlation function 
changes radically as seen in the right plot. Now only every third amplitude 
are correlated and this comes as a result of the period three invariance of the 
model. This plot shows also how the characteristic time-scale changes among 
the different shells. It is seen by the extent of the correlation-peeks in time 
which decreases with shell-number. When relating the characteristic time-scale 
to the turnover-time (r„) this dependence comes direct from dimension analysis 

Hi. 



3.2 Front motion during burst 

The motion of the bursts is a part of a more general motion of different organi- 
sations of the amplitudes travelling with exponential increasing speed from the 
lower towards the higher shells where they vanish because of viscosity [10|. A 
way to see this is to look at the changes in the instantaneous amplitude spectre 
during the motion of a burst. This is shown in Figj^ by snapshots of In |un(i)| 
vs. n where the time between snapshots decrease by a factor of l/\/2 giving 
roughly an equidistant motion of the burst. As for all other bursts Fig.^ reveals 
that the burst travels through the shells as a front keeping the same overall 
shape. Just at the maximum rise of the amplitudes the overall scaling expo- 
nent of the inertial range is a bit lower than the Kolmogorov scaling-law shown 
by the dashed line in Figjsj. Immediately after the last snapshot all the shells 
enters the oscillatory state. 
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Figure 3: Snapshots of log \ un\ during the cascade of a burst 



3.3 Real- valued model 

Due to the invariance in the model the creation and behaviour of the bursts are 
unaffected by the complex phase of the amplitudes. A model in terms of real 
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values will therefore be used in the following analysis: 



+ ^^n^ = -kn{rn+irn+2 - ^Tn-irn+l - ^rn-2r-n-l) + /5„,4 (4) 

having r„ = \un\, no conjugations and "—1" instead of "i" in front of the 



coupling-term |11 



4 Local variables 

From the construction of the model the dynamics of a given shell depends only 
on the instantaneous configuration of the neighbouring shells and it has no 
explicit dependence on the present or past states. If we at first restrict ourself 
to the inertial range neglecting forcing and viscosity the neighbouring shells 
may be seen as a local phase-space of a shell since their configuration through 
the coupling-term exactly determines the instantaneous dynamics (r„) of the 
amplitude r„. To characterise this local phase-space each set of neighbouring 
shells will be called the local shells: Ln = (r„_2, r„_i, r„_|_i, r„+2) of the n^'^ 
shell, and they should not be seen as part of the other amplitudes but rather 
as an isolated set of variables determining 

The configuration of L„ will be described by first choosing the slope of 
\nx{Ln) which is nothing but the local scaling exponent at the n*^ shell. To 
continue we define 

ffn = lnA(L„) (5) 

and choose the mean, curvature and third order component of ffn- This gives 
the local variables: Pn = {An, Bn, Cn, Dn) of r„ defined as the coefficients of 
the projection of ff„ on the orthogonal basis given by the matrix T : 



2V^ 



where 
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(7) 



and a = 1/4 and /? = 1/10. 

The basis of the local variables is plotted in Fig.^ showing how it can be 
characterised as a simple "Taylor-series" expansion of ffn- These variables are 
believed to be the right variables to monitor the dynamics of the model since 
they describe globally the configuration of the local shells instead of focusing on 
the individual neighbouring shells. The local scaling of shell-models has been 
studied earlier |9|, ^ , but this was the instantaneous local scaling averaged over 
all shells and using a coarse-grained time resolution. 



6 



-0 — 


B 


H p- 










\ 

\ 


A 


/ 

/ 


0.2 








\ 




/ 


0.1 


\ 

/ \, 




h- 




\ 


/ 




/ \ 








\ 


/ 

/ 







/ 






\ 

\ 


/ 


-0.1 




\ / 






\ 


/ 






\ / 






\ c 


/ 

/ 


-0.2 


-'*"' B 






















n-2 n-1 n n+1 n + 2 

l,ocal Shells 



n-2 n-1 n n+1 n+2 

Local Shells 



Figure 4: The basis of the local variables 



4.1 Application to the model 

To implement Eq.^J^ into the model we assume the components of fjn to behave 
smoothly in n such that r„ ~ 2"^", giving: 

j. — _ 2^^" _ 5_2 — 'iCn _ 1 — 5 2~3^n+£'n^ — y j^^ 2^" (8) 

Eq.^ gives direct evidence of the period three invariance of the model: Since 
the dynamics only depends on the combinations (3i?„ — C^, A^), we define 
En = 3Bn — Dn- The model is then invariant to the orthogonal component of 
En- -LEn = SDn + Bn which is nothing but the period three invariance as seen 
in in Fig.|5[ 




Figure 5: The period three invariance of the local basis 

From the construction of Eq.^ is should be noted that the sign of r„ and 

thereby the monotony of the dynamics is only a function of En and C„ when 

neglecting the viscosity-term. Because An is outside the brackets it will affect 

th 

the response-time of the dynamics. Now the dynamics of the n amplitude is 
only determined by three local variables: 

= {Em Cn, An). 

Even though this new set of local variables (Vn) form a efficient phase-space it 
should not be confused with the actual 2N-dimensional phase-space of the free 
variables in the model. 
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4.2 The Local Attractor of the model 

Since is a local phase-space the trajectory of Vn{t) in time will describe a 
three dimensional local attractor of the n^^ shell dynamics. Fig.^ show the local 
attractor of the 14*^^ shell during a time-span of two successive bursts where 
some additional features is placed to explain the dynamics of the attractor: 




Figure 6: The local attractor of the 14'th shell and its projection on a {En, Cn)- 
plane together with the surface of r^ = 0, the Kolmogorov fixed-point-line and 
arrows of characteristic flow. 

First we note that the trajectory is projected down on a (E'n, C„)-plane to 
help giving a three dimensional understanding of the attractor. Then we focus 
on the vertical line which correspond to the Kolmogorov fixed-point given by 
{En, Cn, Ajy) = (—1,0, •). Right after every bursts the trajectories encircles this 
line during the relaxations. As the oscillations die out the dynamics slow down 
making the trajectories stay close to the region of ~ in In Fig.^ the 
curved sheet is the manifold of r„ = derived from Eq.^ and it is seen how 
the trajectory stays close to the manifold. ( note that the trajectory is shown 
thinner for negative r„) 

When a burst approaches from the lower shells it affects the configuration of 
local shells forcing the trajectory away from the manifold. This causes rn and 
thereby r„ to increase rapidly making the shell participate the burst. During 
the burst the trajectory approaches the Kolmogorov fixed-point-line around 
which it begins to circle again etc. The same behaviour repeats throughout 
the evolution of the model making the local attractor capture all the general 
dynamics of the model. 
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Every other shell participating in the burst has qualitatively the same local 
attractor with the same characteristics. It should be noted that if the viscous 
term only affects the last shells, abandoning the inertial range, the model would 
still produce bursts and in this case the oscillations would not bend off but follow 
the Kolmogorov fixed-point strait down until the next burst approaches. 

5 The cause of intermittency 

From the behaviour of the local attractor it is possible to explain the inter- 
mittent shift between bursts and oscillatory relaxation of the model. What is 
needed is the answers to the following two questions: Why is the manifold of 
r„ = stable, attracting the oscillatory state into a relaxing period and what 
changes this stability as a burst approaches. 

5.1 Creation of the relaxing period. 

To analyse the stability of the manifold we have to know the flow in the phase- 
space Vn and this will be done by estimating An, En,Cn' 

First we assume again r„ 2^" to get r„ ^ ln(2) 2^"yi„ which will be used 
to estimate A„. Then we insert An into the transformations of EqJ^ getting 
En and C„ as a function of An+j, j = {—2, 1—, 1, 2}. To proceed we note that 
because of the regular dynamics during oscillations all the local variables for 
the different shells are roughly equal despite a Kolmogorov-scaling of the mean 
values (An). This makes us assume the following condition between the local 
variables: 

(An+j, En+j, Cn+j) ~ (^An — i, En, , j = { — 2, 1 — , 1, 2} (9) 

When inserted into the different An+j it causes En and C„ to resemble An 
within pre-factors in front of the coupling- and viscous-terms. 
As result the monotony of En and Cn follows that of An- 

Now the general flow in Vn only depends on the sign of r^i changing at the 
manifold and indicated by the arrows showed in Fig.|6[ From the orientation of 
the flow and the position of the manifold the trajectory is caused to close in on 
the manifold and slowly drift downwards creating a relaxing period. 

5.2 Bursts 

The stability of the manifold and thereby of the relaxing state depends critically 
on the condition of Eq.|9| used in the derivation above. The thing that destroys 
this condition is the approach of a burst from the lower shells, affecting only 
'''n-2,i^n-i- The manifold then loses its stability and the state is forced into a 
region of strong positive r„ making the shell participate in the burst. Now as 
changes violently it causes the manifold of the higher shells to become unstable 
etc. and thus the burst spreads through the shells because of a chain-reaction. 
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6 Conclusion 



In this article the standard GOY shell model has been analysed on the basis of 
its dynamics rather than its statistics. A detailed analysis of the time-evolution 
revels the following: 

The dynamics of the model follows two different states where violent bursts 
are interchanged by an oscillatory relaxing state. It is showed that the dynamics 
of the shells are mutually correlated and the burst travels through the shells 
like a front. Because bursts in the model cascade nearly unaffected through 
the shells in the inertial range, each set of neighbouring shells entering the 
coupling-terms can be seen as local phase-spaces of the corresponding shells, 
and when expressed in a simple "Taylor-series" base their dynamics describe an 
approximative attractor of the model. From the dynamics of the local attractor 
the intermittency of the model is explained. 
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